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A statistically stationary and nearly homogeneous turbulent shear flow is established by 
an additional volume forcing in combination with stress-free boundary conditions in the 
shear direction. Both turbulent energy and enstrophy are stationary to a much better ap- 
proximation than in previous simulations that use remeshing. The temporal fluctuations 
decrease with increasing Reynolds number. Energy spectra and shear-stress cospectra 
show that local isotropy is satisfactorily obeyed at the level of second-order moments. 
However, derivative moments of high-order up to n = 7 yield increasing moments for 
n > 4 for the spanwise vorticity and the transverse derivative of the streamwise velocity 
in the range of Taylor Reynolds numbers 59 < R\ < 99. These findings, which are in 
apparent violation of local isotropy, agree with recent measurements. 



1. Introduction 

Turbulent flows in nature and laboratory are mostly anisotropic on their largest scales. 
A typical situation is one in which an imposed shear gradient sustains the turbulence. 
One fundamental question is the effect of anisotropic large scales on the statistical be- 
haviour on the smaller scales. Kolmogorov (1941) postulated that at sufficiently large 
Reynolds number statistics of small scales become isotropic even if the large scale driv- 
ing is anisotropic. Consequently, cross-correlations, indicated, e.g., by the shear-stress 
cospectrum E xy (k) in the Fourier space, should decay very fast to zero with respect to 
increasing wavenumber. Based on this postulate, Lumley (1967) predicted, for a turbulent 
shear flow with constant mean shear rate S, the rate at which anisotropy vanishes with 
decreasing scale size. If the shear time scale, T s — S^ 1 , is much greater than the viscous 
time scale, t v = (z^/e) 1 / 2 , in other words, if St v <C 1 holds, he showed by dimensional 
arguments that the shear-stress cospectrum follows E xy (k) ~ e 1 / 3 S'fc _7 / 3 . Quantities v 
and e are the kinematic viscosity and the energy dissipation rate, respectively. Note that 
the cospectrum would be exactly zero in a system that is perfectly isotropic on all scales. 
The concept of local isotropy on the level of second order moments seemed to be sup- 
ported by direct numerical simulations (She et al. 1993) and a number of high-Reynolds 
number measurements in turbulent boundary layers (Caughey et al. 1979; Saddoughi & 
Veeravalli 1994). Very recently, however, Kurien & Sreenivasan (2000) found a slower de- 
cay of anisotropies in an atmospheric boundary layer at R\ w 2000. They reported, e.g., 
a scaling of the cospectrum with E xy (k x ) ~ k x 2A . Deviations from the local isotropy 
were also found in a number of new experiments in nearly homogeneous shear flows 
when attention was focused on higher order moments (Garg & Warhaft 1998; Ferchichi 
& Tavoularis 2000; Shen & Warhaft 2000). 

Quantities which are very sensitive to anisotropic contributions at smallest scales are 
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derivative moments of the turbulent velocity components, u\. We denote moments of odd 
orders (k — 1,2,3) by 

s 2k+l (duj dXo ) - {{du , /dXj)2)k+1/2 , (i-i) 

and of even orders (fc = 2, 3) by 

^ ^« / /« \ (( d K/dxj) 2k ) 
K^l^i) = {{d</dX]W • 

The brackets (•) denote the ensemble average and no tensor summations are applied 
for i and j. If the turbulence in the shear flow is locally isotropic odd moments of the 
transverse derivative (i = x,j = y) should decay with R^ 1 - Shen & Warhaft (2000) 
measured a slower decay of the derivative skewness, S^{du' x /dy) <~ R\°' 5 , with a value 
of 0.2 at R\ ~ 1000. The fifth moment S$(du x / dy) remained constant while Sr(du x /dy) 
even increased in the range of Taylor Reynolds numbers between 100 and 1000. These 
results were confirmed in part by similar experiments of Ferchichi & Tavoularis (2000). 
Third order transverse derivative moments of both experiments agreed also with those 
found in homogeneous shear flow simulations at R\ ~ 100 (Pumir 1996). All the results 
indicate that the influence of large scale anisotropics on the statistics of the smallest 
scales does not decay as fast as predicted by dimensional arguments, and in some cases 
may not decay at all. 

This paper is addressed to the numerical investigation of statistical properties in the 
simplest shear flow configuration, the homogeneous shear flow. The flow is characterized 
by the following mean profiles of the velocity components 

(u x )=Sy, ( Uy ) = ( Uz )=0, (1.3) 

where x are streamwise (or downstream), y shear (or wall- normal) , and z spanwise di- 
rections, respectively. The velocity components can be decomposed in a mean fraction 
and a fluctuating turbulent part, Ui = (ui) + u\ for i = x,y and z. Besides finite differ- 
ence methods with shear periodic boundary conditions (Gerz et al. 1989) pseudospectral 
methods in a fully periodic domain using a remeshing (Rogallo 1981; Rogers & Moin 
1987; Pumir 1996; Gualtieri et al. 2000) were mostly used to model a homogeneous shear 
flow. In this method, it was harder to reach a stationary turbulent state. An alternative 
method to overcome this problem was suggested by Schumacher & Eckhardt (2000) very 
recently. The method avoids the remeshing that needs to be applied for the increasingly 
skewed grid of Rogallo (1981) and that corresponds to a time-periodic driving of the 
flow. A statistically stationary state for long times can be maintained by using stress-free 
boundary conditions in y direction and an appropriate body force. Results for third and 
fourth order agreed with those of Rogers & Moin (1987) and of Pumir (1996). Derivative 
skewness was weakly decaying with respect to R\. This paper extends that analysis. We 
systematically investigate the temporal fluctuations of the energy and enstrophy around 
their mean values and compare them with the findings of Pumir (1996) and Gualtieri et 
al. (2000). We show that the level of fluctuations around the mean quantities is much 
smaller and reaches values that are observed in turbulence measurements. Energy as 
well as shear-stress cospectra are studied. Inspired by recent experiments we extend the 
analysis systematically to higher order derivative moments (n < 7) in the range of Tay- 
lor Reynolds numbers accessible. Non-decreasing higher order moments indicate that 
very rare small-scale events may exist. We shall also discuss the problem of statistical 
convergence of higher order moments in detail. 
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Run No. 1 2 3 4 5 



Re s 


500 


1000 


1000 


2000 


2000 


Rx 


59 


79 


87 


95 


99 


L X 1 Ly 


2tt 


4-k 


2tt 


47T 


2tt 


L z /L y 


2tt 


2n 


2tt 


2tt 


2tt 


N x 


256 


256 


256 


256 


256 


Ny 


65 


129 


65 


129 


65 


N z 


256 


128 


256 


128 


256 




2.82 


2.48 


1.62 


1.54 


1.00 


Ay/n 


0.57 


0.50 


0.99 


0.80 


1.60 


A/r? 


17.8 


18.9 


19.8 


20.2 


20.7 


L/V 


36.4 


63.3 


58.4 


80.5 


77.5 



Table 1. List of the presented runs with its parameters: shear Reynolds number Re s , Taylor 
Reynolds number R\, aspect ratios L x /L y and L z /L y , the numbers of grid points N x , N y , and 
N z , the spectral resolution criterion kmaxV, and the ratios of grid spacing in y direction Ay, of 
Taylor microscale A, and of integral length scale L to Kolmogorov length n. 



2. Numerical model 

With length scales measured in units of the gap width L y , and time scales in units of 
S^ 1 , the dimensionless form of the equations for an incompressible Navier-Stokes fluid 
become 

^ + (u-V)u = -Vp+— V 2 u + f , V-u = 0, (2.1) 

where p{x.,i) is the pressure, u(x, t) the velocity field. The shear Reynolds number is 
Re s = SL 2 /i>. In the x and z directions periodic boundary conditions apply. In the 
other direction the flow domain is bounded by two parallel flat surfaces that are assumed 
to be impenetrable and stress-free, 



du x du z 
dy dy 



u v = ~~a~~ = ~H7~ = for y = 0, L y . (2.2) 



As discussed in more detail in a recent paper by Schumacher & Eckhardt (2001), the 
free slip boundary conditions allow for efficient numerical simulations with Fourier modes 
for the velocity components. The equations are integrated by means of a pseudospectral 
technique using a 2/3-rule de- aliasing. Time integration was done by means of a 5-th 
order Runge-Kutta scheme with adaptive time stepping (Hairer et al. 1993). In table 1 
the computational parameters are summarized for all the runs presented. The aspect ratio 
L x /L z and the shear Reynolds number Re s were varied. It is also visible that the spectral 
resolutions in x and y directions are different causing maximum wavenumber components 
of k x , max (= k ZtTnax ) = 2irN x /(3L x ) and k Vtmax = 2-7riV y /3, respectively, after de-aliasing. 
As a criterion for sufficient spectral resolution k max rj > 1 is used (Pope 2000) with 
Kolmogorov length scale r\ = (z/ 3 /e) 1 / 4 and k max = x/8ttN x /(3L x ) for a uniform grid 
case. For the present model with different grid spacing in the three dimensions, it is 
suggested therefore to take k max = {(k 2 x max + k 2 y max + k 2 z max )/Z) 1/2 . Additionally we 
plotted the ratio of grid spacing in shear direction to the Kolmogorov scale, Ay/n, to 
indicate the sufficient spectral resolution. The Taylor microscale A is defined here with 
the downstream root-mean-square velocity 

, (K) 2 ) 1/2 <.rms . . 

((du'jdxyy/ 2 {du'jdx) rms ' 1 
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Figure 1. Mean velocity profiles of the three velocity components for different Reynolds num- 
bers. The inset magnifies the profiles at the boundary y/L y = 1 where deviations from transverse 
homogeneity are present. The thin solid line corresponds with the exact linear profile. 



and the integral length scale L is given by 

1 f°° 

L =7TTT2\ *-«(x + re x )<(x)). (2.4) 

\\ u x) I JO 

The Taylor Reynolds number follows as R\ — u' x rms X/v. The mean shear and turbulence 
are maintained by a suitable body force f(x, t). An almost linear mean profile (u x ){y) = 
(y — 1/2) for y G [0, 1] (in dimensionless form) can be approximated by a finite Fourier 
sum of cosines 

4 ^ cos[(2n + l)7ry] 
{ux){y)- j> 77. — — 7T 2 • (2-5) 

The external forcing f was chosen such that the six modes used in ( |2.5| ) remained constant 
in time, i.e. dRe{u x (q, t)}/dt = for Fourier modes with q = (2n + l)ire y for n = to 5. 
It was shown by Schumacher & Eckhardt (2000) that this forcing, although itself varying 
with respect to time, causes the expected mean profiles of the velocity field components, 
(u y ) = (u z ) =0 and (u x ) — y (see figure [l]). The effects of the free slip surfaces at y = 
and y = L y on the bulk behaviour are much weaker than those of rigid walls, since only 
the wall-normal component is forced to vanish and a nearly homogeneous, statistically 
stationary turbulent state can be established. 



3. Energy balance and power spectra 

Figure |^ shows the temporal evolution of the specific turbulent kinetic energy E(t) 
and the specific enstrophy f2(i) defined as 

E{t) = \{{u' l f) v , m = \{^?)v, (3.1) 
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Figure 2. Temporal behaviour of the specific turbulent kinetic energy E(t) for three different 
values of Taylor Reynolds numbers, R\ — 59 (dotted line), R\ = 87 (dashed line), R\ = 99 
(solid line), (a) Specific turbulent kinetic energy E(t). (b) Specific enstrophy 0(t). 



where (-)y denotes an average over the volume. The vorticity is defined as uj = V x u'. 
The graphs indicate that a statistically stationary state is established in all cases. Ad- 
ditionally we analyzed the ratio of the standard deviations o~(E) and cr(f2) to their cor- 
responding temporal means, (E) t and (£l) t . Both ratios were found to be smaller in 
comparison to the results of remeshing simulations. Pumir (1996) as well as Gualtieri 
et al. (2000) reported ratios of about 50% for both, energy and enstrophy at R\ ~ 90. 
In contrast, we found a decrease of the ratios from about 23% at R\ = 59 to 10% at 
R\ = 99 (see table 2). The large fluctuations may thus be caused by the periodic driving 
due to the remeshing, though this question still remains open. Statistical stationarity was 
further checked by comparing the ratio of turbulent production, P = —(u' x u' y )d(u x )/dy, 
to the energy dissipation rate, e = i , ((du' i /dxj)' 2 ) . The ratio P/e was always found to 
be very close to 1 (see table 2). The shear parameter 5* = Sq 2 /e with q 2 = ((u-) 2 ) re- 
mained nearly unchanged in the five simulations, consistent with a nearly constant ratio 
— (u' x u' y ) I q 2 , which should equal 1/5* in the stationary case. 

Next we investigated the energy dissipation rate as a function of the Taylor Reynolds 
number. The expectation is that the energy dissipation rate at high Reynolds numbers 
is independent of viscosity, and is a constant of order unity when rescaled on the integral 
scale L and the root-mean-square velocity u' x rms , i.e. e = e L/(u' x rms ) 3 . In table 2 the 
results are given for the present data. Quantitiy e seems to saturate at order unity for R\ 
between 80 and 100 which cannot be exceeded with the present grid resolution. A similar 
behaviour was found in the case of homogeneous isotropic turbulence (Sreenivasan 1998). 

As listed also in table 2, about half of the total amount of velocity fluctuations is 
contained in the downstream component u' x . The downstream fluctuations ((u' x ) 2 )(t) were 
found to have always a similar temporal shape as the total kinetic energies while the mean 
square moments for the both remaining components u' y and u' z have smaller variations 
with respect to time. We observed also that the fraction ((u x ) 2 )/q 2 with q 2 = ((u^) 2 ) 
decreased slowly with respect to R\. Both aspects seem to be connected to the stronger 
fragmentation of coherent streaks and vortices (Schumacher & Eckhardt 2000, 2001) that 
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Run No. 


1 


2 


3 


4 


5 


<«) 2 >/a 2 
{{<) 2 )/ q 2 

~(u' x u' y )/q 2 


0.57 


0.52 


0.55 


0.51 


0.51 


0.12 


0.16 


0.15 


0.18 


0.18 


0.31 


0.32 


0.30 


0.31 


0.31 


0.11 


0.12 


0.11 


0.12 


0.12 


S* 


7.32 


8.18 


8.58 


8.22 


8.37 


St v 


0.38 


0.25 


0.25 


0.19 


0.19 




1.04 


1.00 


0.99 


1.01 


1.01 


e 


1.18 


0.87 


0.70 


0.77 


0.72 


°{E)/(E) t 


0.20 


0.12 


0.11 


0.09 


0.10 




0.23 


0.13 


0.12 


0.09 


0.10 



Table 2. Physical properties of the homogeneous shear flow runs. Compared are the Reynolds 
stress components, {u^Uj), shear parameter, S* , the ratio of Kolomogorov time t v to shear 
time S , the ratio of turbulence production to dissipation, P/e and the dimensionless energy 
dissipation rate e. The last two rows list the ratio of standard deviations to the temporal means 
for the specific turbulent kinetic energy, E — E/(SL y ) 2 , and the enstrophy, f2 = Q./S 2 . 
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Figure 3. Normalized energy spectra and shear-stress cospectra for different values of R\. (a) 
Energy spectra E(k)/(eu 5 ) 1 ^ 4 ' are plotted versus kr/. Straight solid line indicates scaling with 
k~ 5 ^ 3 and Ck = 2. (b) Shear-stress cospectra —E xy (kx)/(e 3 S~ 5 ) 1/l2 are shown versus k x Ls with 
L s = (eS 



-3U/2 



The straight solid line indicates scaling with kx 7 ^ 3 and C xy = 0.15. Argument 



k x was taken here in order to compare recent data to experimental findings. 



is observed with increasing Re and the reduced downstream correlation of the velocity 
fluctuations. 

In figure ^a we plot the kinetic energy spectra for five different Taylor Reynolds num- 
bers. With increasing values of R\ a Kolmogorov scaling range can be observed for inter- 
mediate wavenumbers. The straight line is fitted with a (Kolmogorov) constant Ck = 2. 
It exceeds slightly the values of Ck — 1.5 - 1.6 known from high- Reynolds number mea- 
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Run No. 



S s (lu z ) -7.9 x 1CT 1 -7.1 x KT 1 -6.7 x KT 1 -6.9 x KT 1 -6.3 x KT 1 

K 4 (lu z ) 5.6 x 10° 6.4 x 10° 6.6 x 10° 6.6 x 10° 7.3 x 10° 

S 5 (w») -8.9x10° -1.7 x 10 1 -1.7 x 10 1 -1.8 x 10 1 -1.9 x 10 1 

K 6 (u x ) 4.9 x 10 1 1.3 x 10 2 1.5 x 10 2 1.5 x 10 2 1.9 x 10 2 

S 7 (w*) -2.2 x 10 2 -7.4 x 10 2 -8.6 x 10 2 -8.8 x 10 2 -1.1 x 10 3 



S 3 (du' x /dy) 9.6 x KT 1 9.2 x 10" 1 9.0 x 10" 1 9.2 x 10" 1 8.7 x 10 

K 4 (du' x /dy) 5.6 x 10° 6.4 x 10° 6.4 x 10° 7.2 x 10° 7.1 x 10 

S 6 (du' x /dy) 1.0 x 10 1 2.0 x 10 1 2.0 x 10 1 2.5 x 10 1 2.1 x 10 1 

K 6 (du x /dy) 5.0 x 10 1 1.3 x 10 2 1.3 x 10 2 2.0 x 10 2 1.7 x 10 2 

S 7 (du' x /dy) 2.1 x 10 2 7.4 x 10 2 7.5 x 10 2 1.4 x 10 3 9.5 x 10 2 

S 3 (du x /dx) -7.4 x 10 _1 -4.4 x 10 _1 -5.0 x 10 _1 -2.7 x 10" 

K 4 (du x /dx) 6.3 x 10° 4.8 x 10° 5.3 x 10° 4.3 x 10° 

S 5 (du x /dx) -1.6 x 10 1 -6.9 x 10° -9.2 x 10° -3.7 x 10' 

K 6 (du' x /dx) 1.2 x 10 2 5.8 x 10 1 7.8 x 10 1 4.2 x 10 1 6.8 x 10 1 

S 7 (du' x /dx) -6.1 x 10 2 -1.5 x 10 2 -2.6 x 10 2 -6.9 x 10 1 -1.7 x 10 2 

Table 3. Derivative moments of orders n = 3, 4, 5, 6 and 7 of the transverse derivative du' x /dy, 
the longitudinal derivative du' x /dx and the spanwise vorticity uj z . 



surements (Pope 2000) which can be related to the bottleneck phenomenon (Falkovich 
1994). This effect causes an energy pileup that can result in a larger Ck for the small 
scaling ranges observed here. In figure |^b the cospectra are compared. It is observed that 
their magnitude is always smaller by about two orders of magnitude than the kinetic en- 
ergy spectrum over the whole range of wavenumbers, except the smallest wavenumbers. 
Both shear-stress cospectra fit to the power law proposed by Lumley (1967) in a small 
range of wavenumbers. The constant C xy was found here to be C xy ~ 0.15, equal to the 
values found in high- Reynolds-number experiments (Saddoughi & Veeravalli 1994). As 



already mentioned, recent high-Reynolds number measurements Kurien & Sreenivasan 
(2000) indicated deviations from the k x ' 1 ^ decay law. We note here that our Taylor 



Reynolds number is too small to draw any robust conclusion about the strength of de- 
viations from the classical Kolmogorov-Lumley scaling, but the results do not seem to 
contradict the concept of local isotropy on the level of second order moments. 



4. Higher order derivative moments 

In this section, the results for higher order moments are presented. The quantities 
studied are those derivative moments that display sensitivity to deviations from local 
isotropy. Three different quantites were investigated: the transverse derivative of the 
turbulent streamwise velocity, du' x /dy, the spanwise vorticity, lu z = du' y /dx — du' x /dy, 
and the longitudinal derivative of the turbulent streamwise velocity, du' x /dx. In table 3 
we summarize the results for orders 3 to 7 for five different simulations. The moments 
of the longitudinal derivative du' x /dx do not vary in order of magnitude for all values 
of n, and are found to agree with Ferchichi & Tavoularis (2000) as well as with Shen & 
Warhaft (2000). For a larger range of Taylor Reynolds numbers both experiments find a 
slow increase for all longitudinal moments. If the turbulent flow would be isotropic, the 
odd moments S n (du' x /dy) would be exactly zero for all odd n > 0. In contrast, we observe 
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Figure 4. Higher order moments of the transverse derivative du' x /dy are compared with other 
numerical and experimental findings, (a) Skewness S3, (b) Flatness K4. (c) Hyperskewness S5, 
hyperflatness Kq and the 7th moment SV are compared. □, Pumir (1996); A, Shen & Warhaft 
(2000); present data. Error bars for are calculated with respect to the variation of the data 
points in the y direction. 
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Figure 5. Statistical convergence of higher order derivative moments. Left column: Results for 
the normalized transverse derivative are shown, probability density function p(z), z p(z) and 
z 7 p(z). Right column: The same functions for the normalized longitudinal derivative are plotted. 
Both data sets are for R\ = 95 (run 4). 
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nonzero values for the transverse derivative and spanwise vorticity moments. While the 
third order slowly decreases in the range of Taylor Reynolds numbers, fifth and seventh 
order increase with R\ . In figure ^ we compare our results with the numerical data of 
Pumir (1996) for Sz{du' x /dy) and with the experimental data of Shen & Warhaft (2000) 
at R\ ~ 100. The error bars for our data display the standard deviation of the y profiles 
of the moments. Averaging was always performed in time and in planes at fixed y. Data 
points closest to the boundaries are excluded from the evaluation of the averages. The 
snapshots of the full velocity field were always separated by at least one shear time unit, 
St = 1. The minimal number of used snapshots was 114 (for run 5). The results are in 
agreement with previous numerical and experimental findings. 

The fluctuations around the mean grow for increasing order as expected. These fluc- 
tuations were subject to a more detailed investigation. Figure ^| shows the probability 
density functions (pdf), denoted asp(z), of z — (du' x /dy)/((du' x /dy) 2 ) 1 ^ 2 in the left col- 
umn and of z = (du' x / dx) / ((du' x / dx) 2 ) 1 / 2 in the right column. The results are calculated 
from 6 x 10 8 data points at R\ = 95. A check of the statistical convergence of the nth 
order moments can be made by plotting the integrands z n p(z) over z which is also shown 
in figure |[ While the longitudinal derivative has converged well, the integrands for the 
transverse derivative show a strong scatter in the tails of the pdf. The area that covers 
the scattered data measures the uncertainty around the mean in an alternative way. 

In the following, we wish to discuss the strong scattering of higher order moments 
and will demonstrate by means of a simple analytical model that this is a more general 
feature for data analysis. The pdf (see figure ||) can be roughly approximated by the 
following normalized exponential distribution p(z) = (acosh(o;2:o)/2) exp(— a\z — z Q \) for 
which all moments can be evaluated analytically, following to (z n ) ~ nl . Let N be the 
number of statistically independent samples. The standard deviation, which is defined as 
a{z n ) = [{(z 2n ) - (z n ) 2 )/(N - l)] 1 / 2 , follows then as 



(2tQI - (tanh(c^ )n!) 2 l 1/2 

N ' [ ' ' 



for the odd normalized moments S n = J z n p(z)dz. The denominator was simplified 
because N ^> 1. When taking, e.g., a = 1 and zq = 1, the standard deviation for 
the seventh moment decreases from 24% at N = 10 5 to 1% at N = 10 s . Thus, the 
calculation demonstrates clearly the sensitivity of higher order moments on the number 
N which can give rise to very large fluctuations around the mean as observed in our 
simulations. While moments itself are proportional to the factorial of the order, nl, its 
fluctuations become even proportional to (^n)!) 1 / 2 . Tennekes & Lumley (1972) estimated 
N(l) = T/(2r) for a first order moment, where r is the integral time scale and T the 
duration of the measurement. Starting from this point, Sreenivasan et al. (1978) extracted 
relations for the accuracy of higher order moments using their atmospheric boundary 
layer data. Following their lines, the number of statistically independent samples for an 
nth moment, N(n), would result in N(n)/N(l) = (0.82 - O.Q7n) -1 . Aside from the fact 
that the fluctuations for fixed N and growing n can increase rapidly, this dependence 
would improve the statistical convergence weakly. 

In summary, we have presented an alternative method of modelling a statistically sta- 
tionary homogeneous shear flow by using stress-free boundary conditions in y direction. 
On the level of second order moments our results are not in conflict with the classical 
concepts of local isotropy as indicated by the power spectra of the velocity correlations. 
We found non-decreasing derivative moments for the orders 4 to 7 which agree with the 
experimental findings of Shen & Warhaft for R\ ^100 but are in apparent violation 
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of local isotropy. A next step would be to relate typical structures of turbulent shear 
flows, such as streamwise vortices and streamwise streaks, to the statistical properties 
presented here. 

The author thanks B. Eckhardt, A. Pumir, K. R. Sreenivasan, P. K. Yeung and 
Z. Warhaft for fruitful discussions and helpful comments. The numerical simulations 
were done on a Cray T-90 at the John von Neumann-Institut fur Computing at the 
Forschungszentrum Jiilich. This work was sponsored in parts by the Feodor-Lynen Fel- 
lowship Program of the Alexander von Humboldt -Foundation and Yale University. 
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